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The macroscopic behavior of many physical systems can be approximately described by classical 
quantities. However, quantum theory demands the existence of omnipresent quantum fluctuations 
on top of this classical background - which, albeit small, should have some impact onto its dynamics. 
The correct treatment of this quantum back-reaction is one of the main problems in quantum gravity 
and related to fundamental questions such as the initial (big bang) singularity or the cosmological 
constant. By means of the qualitative analogy between gravity and fluid dynamics, we try to shed 
some light onto these problems and show some of the difficulties associated with the calculation of 
the quantum back-reaction starting from the classical (macroscopic) equation of motion. 



In physics (and other sciences), it is often useful to make gedanken experiments in order to explore the limiting 
cases of a given theory and to approach it from different sides. Therefore, let us do the same here and imagine the 
following situation: While sitting on the coffee table in your physics department, a colleague approaches you to ask 
you a favor: "Hi, I have a question for you - there is this theory with the following properties:" 

• The classical equations of motion are known and my experimental collaborators have them reasonably well 
tested at large length scales and low energies. 

• However, we all agree that this classical description should break down at small length scales, where somehow 
quantum effects should become important. 

• But if I try to quantize this theory in the usual perturbative manner, I get non-renormalizable UV-divergences 
and apparently I have to introduce new (unknown) parameters. 

• Unfortunately, I do not know and understand the the full quantum theory (yet). 

"Nevertheless, my experimental collaborators are now preparing precision measurements and they want me to estimate 
the back-reaction of the quantum fluctuations onto the classical background profile (as described by the classical 
equations of motion) . Do you have any idea how to do this?" 



From reading the above points, one might have expected that this colleague works in quantum gravity - but in 
fact, precisely the same problem does also apply to fluids dynamics: 

• The classical description is provided by the Euler equation 



• However, this equation clearly breaks down at small length scales, where the continuum fluid picture fails. 



MOTIVATION 



FLUID DYNAMICS 




(1) 



together with the equation of continuity p + V • (pv) — plus the equation of state p = p{p) and the external 
force density / ext . 
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• The perturbative quantization of the longitudinal degrees of freedom yields phonons as quasi-particles and 
the inclusion of phonon loops produces non-renormalizable UV-divergences. Moreover, the quantization of the 
transversal degrees of freedom, i.e., the vorticity Vx», cannot be inferred from the Euler equation at all. E.g., 
in superfluids, one can only have an integer number of vortices and the circulation quantum of a single vortex 
depends on the mass m of the constituent particles - which is not contained in Euler equation. 

• Even in those situations where the the full quantum theory (i.e., the many-particle Hamiltonian) is in principle 
known, it is practically impossible to solve - except in a few simple cases. (We shall study one of them below, 
see also [l[.) 

FIRST TRY 

Clearly, the phonon modes derived from Eq. (fTJ) must also possess quantum fluctuations and these quantum fluctu- 
ations should somehow react back onto the fluid. Let us try to estimate this impact. In our first attempt, we assume 
zero knowledge of the microscopic structure of the fluid and just use the classical equation of motion. In order to avoid 
the aforementioned problems with vorticity, we assume an irrotational flow as described by the Bernoulli equation 

h^ ( f > f + v^ t = -h{Q) = -g Q , (2) 

with <fi being the velocity potential v — V0. The specific enthalpy dh = dp/ g determines the equation of state, where 
we have chosen the relation h(g) = gg describing Bose-Einstein condensates for later comparison. 

Now the usual procedure consists in replacing the c-number fields p and (f> in the equation above by operators 
(more precisely, operator-valued distributions) and to interpret their expectation values as the classical contribution 
(mean-field expansion) 

4> -► <j> = (<j>) + 5$ = (f> d + 54> 

Q -> g= (g) +dg= g cl + 5g. (3) 
Inserting the above split into the Bernoulli equation and taking the expectation value yields 

0d + i(W cl ) 2 + Kxt + <7&l = ™((V<ty) 2 ) , (4) 

where the last term can be attributed to the quantum fluctuations of the phonon modes and should include their 
zero-point pressure. The quantum fluctuations 6(j> of sound derived from Eq. (|2|) behave exactly like a scalar (Klein- 
Fock-Gordon) quantum field - provided that one replaces the speed of light with the sound velocity c s in all expressions. 
Hence the expectation value in the above equation is UV-divergent and behaves as J kd 3 k. Remembering that the 
fluid picture is only valid for large length scales, one would cut off the ^-integration at some maximum wavenumbcr 
fc cut . Below this cut-off k cu t, our description should work fairly well and hence one would arrive at the following 
estimate (see also [H) 

Pzcro-point = -^((V<50) 2 ) OC k* ut ■ (5) 

An analogous line of reasoning has been adopted for estimates of the cosmological constant (i.e., the gravitational 
impact of the quantum vacuum) A tx Mp ianck , which turns out to be many orders of magnitude too large. However, 
in the above approach the main contribution to the zero-point pressure stems from phonon modes with wavenumbers 
close to k cu ti where the theory breaks down. Therefore, one might worry whether the above estimate is reliable. 

SECOND TRY 

In order to tackle this question, let us make a second attempt by incorporating some knowledge about what is 
going on at high wavenumbers k cu t, i-e., how the fluid picture changes at small length scales. Of course, this behavior 
will depend on the kind of fluid under consideration - i.e., we have to specify it. In the following, we shall consider 
dilute Bose-Einstein condensates, since we understand the microscopic structure of these quantum fluids quite well. 
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In Bose-Einstein condensates, the Bernoulli equation breaks down - i.e., acquires additional terms - at the healing 
length £ = l/^/rngg (h = 1 throughout) 

<fi + x(V0) 2 + Vext + 9Q — 77— . (6) 

2 2m ^fg 

For large length scales k <C l/£, the last termQ can be neglected and we recover the usual sound waves with a linear 
dispersion relation uj = c s k. For higher wavenumbers, however, the dispersion relation changes and approaches the 
energy-momentum relation ui — k /(2m) for free particles (with mass m) in the limit kt; 3> 1 

If we now insert the same mean-field expansion §3§ as in the previous Section into the modified Bernoulli equation 
and neglect terms of third or higher order in Sg, we get 

Pzcro-point = ™<(^ Hf) + ^ <(W£) 2 > + 0(Sg 3 ) . (8) 

Without respecting the altered dispersion relation |(7J), the additional term ((WSg) 2 ) would be even stronger UV- 
divergent J k 3 d 3 k than the first one ((V<50) 2 ). However, including the full dispersion relation |(7J), both expectation 
values behave as / k 2 d 3 k for large k and the leading UV-divergences cancel exactly. Nevertheless, the sub-leading 
UV-singularities remain and one obtains 

Pzexo— point gk cu ^ . (9) 

Therefore, the k^ nt behavior in estimate © we just "successfully" calculated in the previous Section is completely 
canceled by the additional term in Eq. ([5]), which becomes relevant for wavenumbers k — 0(l/£) where the fluid 
picture breaks down. Furthermore, even after fully incorporating the cut-off fc^ut ~ l/£ stemming from the healing 
length, the result is still UV-divergent (though the singularity is weaker now). Naturally, one would expect that this 
remaining divergence k 3 ut reflects the omitted details of the microscopic structure of two-particle interaction potential 
which generates the internal pressure in the condensate, cf. 

THIRD TRY - FULL THEORY 

Fortunately, for dilute Bose-Einstein condensates, we are able to address this question analytically. The full many- 
particle Hamiltonian describing the condensate reads 



H = d 3 r 



+ d 3 rd 3 r'i>\r)&(r')V iu t(r-r')y(r)4>(r'), (10) 



(wt). ( w) ■ 

2m 

where Vi n t(r — r') denotes the two-particle interaction potential. In the s-wave approximation, the microscopic 
structure of this potential is omitted and the Hamiltonian simplifies to 



H = d 3 r 



2m 2 



(11) 



In order to make contact to the results of the previous Sections, we introduce the density g — and phase <j> 

operators via the quantum Madelung ansatz ^ = e 1 ^ \fg. Note, however, that this split is highly singular due to the 
commutator \^>(r), W (f ')] = S 3 (r — f '). Therefore, if we try to re-express the above Hamiltonian in terms of g and <j> 

(#t)2|,2 = ($t$)2 _ ^t^^rtjf = g 2 - g6 3 (0) , (12) 

we encounter UV-divergent operator-ordering remnants 5 3 (0) oc k 3 ut . Thus the quantum analogue of the Bernoulli 
equation 

^=i[H,4>] = ^-V ext -gg+^5 3 (Q), (13) 
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contains at least one singular term g6 3 (0) oc gk 3 ut . Comparison with Eq. ([9]) shows that the estimate of p Z oro-point 
in the previous Section is completely spoilt by such unphysical operator-ordering remnants. One might argue that a 
constant (though infinite) pressure shift does not play any role - but the coupling constant g may well depend on an 
external magnetic field and thereby vary spatially (which would feign a pressure gradient). 

Note that the singularity g5 3 (0) oc gk 3 ut could be cured by a smooth two-particle interaction potential Vi nt (r — r '). 
However, insertion of the quantum Madelung split i& — e 1 ^ ^/g into the kinetic term (V^t) . (V^) again yields 
operator-ordering singularities which are independent of Vi nt (r — r') and thus the quantum Bernoulli equation above 
is still ill-defined. Therefore, in order to calculate the back-reaction of the quantum fluctuations correctly, one should 
not use the ill-defined phase operator <p but employ the fundamental field operators i&* and $ instead. In analogy to 
Q , the mean-field split W = ip c + \ into the condensate order parameter tp c plus small quantum fluctuations (phonon 
modes) \ yields the Gross-Pitaevskii equation 

V>c + .g(x 2 )C (14) 

including the quantum back-reaction terms (x'x) and (x 2 )- The first term (x'x) i s called the quantum depletion and 
is finite even in the s-wave approximation Vi n t(f— r') — > gS 3 (f — r '). The second term (x 2 ) (known as the anomalous 
contribution) would be UV-divergent in this approximation Vi n t(r — f ') — > g5 3 (r — f '), but only weakly (x 2 ) ^cut- 
Taking into account the finite range of the interaction Vi n t(r — f ') ^ gS 3 (r — r '), however, it is also finite. Therefore, 
the description in terms of the fundamental field operators W and ^> is well-defined and finite - whereas the phase 
operator <f> is ill-defined and generates unphysical singularities. 



iipc 



2m 



+ V c ^+g\^ c \ 2 + 2g(x f x) 



CONCLUSIONS 



The presented gedanken experiment is based on a qualitative analogy between (quantum) fluids and (quantum) 
gravity sketched in the following table. Even though the problem of quantizing gravity is roughly as old as the quest 
for a full quantum description of fluids, our progress in the latter case is far more explicit (see the various question 
marks below). This difference can probably be mostly attributed to the vastly different input from the experimental 
side. 

Classical equation of motion: Euler equation (JTJ) Einstein equations 

Quantized linearizations: Phonons x Gravitons (?) 

First UV cut-off scale: Healing length £ «-> Planck scale Mpi anc i;? 

Full quantum theory: Many-body Hamiltonian (|10|) «-> Quantum gravity? 

Quantum back-reaction: Zero-point pressure p Z ero-point Cosmological constant A? 

Now, what can we learn from our gedanken experiment, i.e., what are possible lessons for quantum gravity (based on 
our superior understanding of quantum fluids)? 

• Usually, quasi-particles are obtained by quantizing the linearized fluctuations of the classical equation of motion. 
For phonons, this procedure indeed yields the correct description for low energies (on the linear level). However, 
the vorticity cannot be quantized in this way and requires some knowledge of the microscopic structure of the 
fluid - which is not contained in the Euler equation. Whether this linearized quantization procedure works for 
gravitons is not completely clear. One might expect that they behave like Goldstone modes and thus are similar 
to phonons - but without knowing the full quantum theory, this is not evident. 

• Going beyond linear order, one encounters (non-renormalizable) UV divergences in both cases. As we have 
observed in the previous Sections, naively summing the zero-point fluctuations of the phonon modes up to 
some cut-off (and thereby extrapolating the Euler equation to large k) does not give the correct result. This 
observation provides a different view on the naive estimate A oc Afp lanck of the cosmological constant, cf. [J]. 

• A third point is the existence of two (or more) vastly different UV cut-off scales. In Bose-Einstein condensates, 
the first deviation from the Bernoulli equation occurs for length scales at the healing length fc^ ut ~ l/£, where 
the dispersion relation changes. In gravity, a change of the dispersion relation corresponds to the breakdown or 
modification of Lorentz invariance. However, even including the modified dispersion relation ([7]), there are still 
UV divergences in the s-wave scattering approximation pip , which can be removed by taking into account the 
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full two-particle interaction potential VS nt (r — ?'). Therefore, the range of Vj llt (r — r') introduces another UV 
cut-off fc^ ut which is (in dilute Bose- Einstein condensates) much larger fc^ ut 3> fc^ut ~ l/£ than the wavenumbers 
at which the effective Lorentz invariance breaks down. This observation poses the question of the significance 
of the UV scales in quantum gravity - does really everything happen at the Planck scale and is there indeed 
nothing going on at higher scales? 

• In dilute Bose-Einstein condensates with a pure contact interaction, the first cut-off fc^ut ~ l/£ = ^/mgg is 
directly related to the circulation quantum 

r= Idf-v = — N = 2tt£c s N, (15) 
/ m 

which determines the "quantization" of vorticity. In other superfluids, however, the two scales can be indepen- 
dent. In the presence of finite-range (e.g., dipolar) interactions, for example, the first deviation of the dispersion 
relation is determined by the maxon hill/roton dip, which depends on the interaction structure - while the 
circulation quantum is still determined by the mass m. (I.e., in this case, there are three UV scales.) 

• Finally, the presented calculations clearly demonstrate the importance of identifying the correct fundamental 
variables. Given a classical equation valid at large scales, one can generate almost every outcome for the impact 
of the quantum fluctuations by quantizing it in terms of different variables and using the operator ordering 
ambiguities. In fact, the power of fluid dynamics partly relies on the fact that many systems which are totally 
different on the microscopic level obey very similar classical equation of motions at large length scales and low 
energies. Therefore, it is crucial to identify universal (emergent) features and to distinguish them from system 
specific properties. 
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